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Abstract. - We study experimentally the fluctuations of the injected and dissipated energy in 
a system of a colloidal particle trapped in a double well potential periodically modulated by 
an external perturbation. The work done by the external force and the dissipated energy are 
measured close to the stochastic resonance where the injected power is maximum. We show that 
the steady state fluctuation theorem holds in this system. 
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A colloidal particle, confined in a double well poten- 
tial, hops between the two wells at a rate r^, named 
the Kramers' rate, which is determined by the height SU 
of the energy barrier between the two wells, specifically 
r'fc ~ T^^ cxp{~-^Hj;), where Tq is a characteristic time, 
ks the Boltzmann constant and T the heat bath temper- 
ature [1] . When the double well potential U is modulated 
by an external periodic perturbation whose frequency is 
close to rfc the system presents the stochastic resonance 
phenomenon [2], i.e. the hops of the particle between 
the two wells synchronize with the external forcing. The 
stochastic resonance has been widely studied in many dif- 
ferent systems and it has been shown to be a bona fide res- 
qnance looking at the resident time [2,3] and the Fourier 
transform of the signal for different noise intensity [4] . Nu- 
merically, the stochastic resonance has been characterized 
by computing the injected work done by the external agent 
as a function of noise and frequency [5,6]. However the 
fluctuations of the injected and dissipated power at the 
stochastic resonance have never been studied experimen- 
tally. This is a very important and general issue within the 
context of Fluctuation Theorems (FT) which constitute 
extremely useful relations for characterizing the probabili- 
ties of observing entropy production or consumption in out 
of equilibrium systems. These relations were first observed 
in the simulations of a sheared fiuid [7] and later proven 
both for chaotic dynamical systems [8] and for stochastic 
dynamics [9]. These works lead to different formulations 
which find powerful applications for measuring free-energy 



difference in biology (sec e.g. [10] for a review). The hy- 
pothesis and the extensions of fiuctuation theorems [11] 
have been tested in various experimental systems such 
as colloidal particles [12-14], mechanical oscillators [15], 
electric circuits [16] and optically driven single two-level 
systems [17]. The effect of nonharmonic potential on the 
motion of a colloidal particle has been tested by Blickle et 
al. [13,17,18]. However the kind of nonharmonic poten- 
tial used in these experiments did not induce a bistable 
dynamics of the Brownian particle, and as far as we know 
there is only one experimental study of FT for bistable 
systems which are not thermally activated [17]. Recently 
a numerical study, which has explored the distributions 
of the dissipated heat and of the work in a Langcvin dy- 
namics near the stochastic resonance, has shown that FT 
holds in the long-time limit [19,20]. 

To give more insight into this problem we study exper- 
imentally the Steady State Fluctuation Theorem (SSFT) 
in the case of a colloidal bistable system, composed by a 
Brownian particle trapped in a double well potential peri- 
odically modulated by an external driving force. We mea- 
sure the energy injected into the system by the sinusoidal 
perturbation and we analyze the distributions of work and 
heat fluctuations. We find that although the dynamics of 
the system is strongly non-linear the SSFT holds for the 
work integrated on time intervals which are only a few 
periods of the driving force. 

The experimental setup is composed by a custom built 
vertical optical tweezers made of an oil-immersion objec- 
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tive (63x, N.A.=1.3) which focuses a laser beam (wave- 
length A = 1064 nm) to the diffraction limit for trap- 
ping glass beads (2 fiia in diameter). The silica beads 
are dispersed in bidistilled water in very small concentra- 
tion. The suspension is introduced in the sample chamber 
of dimensions 0.25 x 10 x 10 mm^, then a single bead is 
trapped and moved away from others. The position of the 
bead is tracked using a fast-camera with the resolution 
of 108 nm/pixels which gives after treatment the position 
of the bead with an accuracy better than 20 nm. The 
trajectories of the bead are sampled at 50 Hz. 

The position of the trap can be easily displaced on the 
focal plane of the objective by deflecting the laser beam 
using an acousto-optic deflector (AOD). To construct the 
double well potential the laser is focused alternatively at 
two different positions at a rate of 5kHz. The residence 
times Ti (with i = 1,2) of the laser in each of the two po- 
sitions determine the mean trapping strength felt by the 
trapped particle. Indeed if ti = T2 = lOOfis the typi- 
cal diffusion length of the bead during this period is only 
5 nm. As a consequence the bead feels an average double- 
well potential: Uo{x) = ax^ — bx^ — dx, where a, b and d 
are determined by the laser intensity and by the distance 
of the two focal points. In our experiment the distance 
between the two spots is 1.45 /im, which produces a trap 
whose minima are at Xmin = ±0.45/^r?T,. The total inten- 
sity of the laser is 58 mW on the focal plane which cor- 
responds to an inter-well barrier energy SUo = 1.8 ksT. 
Starting from the static symmetric double-trap, (ti ~ T2) 
we modulate the depth of the wells at low frequency by 
modulating the residence times (r,) during which the spot 
remains in each position^ . The modulation of the average 
intensity is harmonic at frequency / and its amplitude 
(t2 — ti)/{t2 + Ti), is 0.7 % of the average intensity in the 
static symmetric case. Thus the potential felt by the bead 
has the following profile in the axial direction: 

U{x,t) = UQ{x) + Up{x,t) = UQ + cx sin(27r/t), (1) 

with aa;^jj„ = 1.8 ksT, bx'^^^ = 3.6 fcsT, d\xmin\ = 
0.44 ksT and c|a;mm| = 0.81 ksT. The amplitude of the 
time dependent perturbation is synchronously acquired 
with the bead trajectory.^ 

An example of the measured potential for t ^ jj and 
^ is shown on the Fig. la). This figures is obtained by 
measuring the probability distribution function P{x, t) of 
X for fixed values of csin(27r/f), it follows that U{x,t) = 
-ln(F(x,t)). 

The X position of the particle can be described by a 
Langevin equation: 



jx ■ 



dU 
dx 



-e, 



(2) 



^We keep the total intensity of the laser constant in order to 
produce a more stable potential. 

^The parameters given here are average parameters since the co- 
efficients a, b and c , obtained from fitted steady distributions at 
given phases, vary with the phase (Sa/a Ri fO%, Sb/b Ri Sc/c Ri 5%). 



with 7 = 1.61 10 ^N s m ^ the friction coefficient and ^ 
the stochastic force. The natural Kramers' rate (c = 0) for 
the particle is Vk = 0.3Hz at T = 300K. When c 7^ the 
particle can experience a stochastic resonance when the 
forcing frequency is close to the Kramers' rate. An exam- 
ple of the sinusoidal force with the corresponding position 
are shown on the figure lb). Since the synchronization is 
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Fig. 1: a) The perturbed potential at f = i and half a forcing 
period later, b) Example of trajectory of the glass bead and 
the corresponding perturbation at / = 0.1 Hz. 



not perfect, sometimes the particle receives energy from 
the perturbation, sometimes the bead moves against the 
perturbation leading to a negative work on the system. 

In the following, all energies are normalized by fcsT. 
From the trajectories, we compute the stochastic Wg and 
the classical Wd works done by the perturbation on the 
system and the heat Q exchanged with the bath. These 
three quantities are defined by the following equations as 
inref. [21]: 



Ws[x{t)] 



Wci[x{t)] 



to 



*"+*^ dU{x,t) 

'^"\tx'J^ 
ox 



(3) 



j-ta+tf 

Q[x{t)] = dx 

Jtn 



dx 



where in this case i/ = j is a multiple of the forcing pe- 
riod. We use both Wg and Wd because they give com- 
plementary informations on the fluctuations of the en- 
ergy injected by the external perturbation into the system 
(see ref. [22] and reference therein for a discussion on this 
point). For example, as we will see later, Ws/T is the to- 
tal entropy production rate in this specific case [23] . The 
heat and the work, defined in eq.3, are related through the 
first principle of thermodynamics: Q = —AU + Ws, where 
At/ — U{x{tf+to),to + tf) — U{x{to),to), whereas the two 
works are related by a boundary term W = —AUp + Wg, 
where AUp = Up{x{tf + to),*/ + ^o) - Up{x{to),to). We 
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identical fluctuations, which will be studied in detail in 
the following of the article. We now focus on the distri- 
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Fig. 2: Injected energy in the system over a single period as 
a function of the driving frequency {Wa □ and Wd o). The 
error bars are computed from the standard deviation of the 
mean over different runs. Inset: Standard deviations of work 
distributions over a single period normalized by the average 
work as a function of the frequency (same symbols) . 



first measure the average work received over one period for 
different frequencies (i/ = i in eq. 3). Each trajectory is 
here recorded during 3200 s in different consecutive runs, 
which corresponds to 160 up to 6400 forcing periods, for 
the range of frequencies explored. In order to increase the 
statistics we consider 10'"' different to- The figure 2 shows 
the evolution of the mean work per period for both defini- 
tions of the work. First, the input average work decreases 
to zero when the frequency tends to zero. Indeed, the bead 
hops randomly several times between the two wells during 
the period. Second, in the limit of high frequencies, the 
particle has not the time to jump on the other side of the 
trap but rather stays in the same well during the period, 
thus the input energy is again decreasing when increasing 
frequency. In the intermediate regime, the particle can 
almost synchronize with the periodical force and follows 
the evolution of the potential. The maximum of injected 
work is found around the frequency / « 0.1 Hz, which is 
comparable with half of the Kramers' rate of the fixed po- 
tential tk = 0.3 Hz. This maximum of transferred energy 
shows that the stochastic resonance for a Brownian par- 
ticle is a bona fide resonance, as it was previously shown 
in experiments using resident time distributions [3,24] or 
directly in simulations [5,6]. It is worth noting that the 
average values of work in this case do not depend on their 
definitions: only the boundary terms, which vanish in av- 
erage with time, are different. In the inset of figure 2, 
we plot the normalized standard deviation of work distri- 
butions [a I lyV) ) as a function of the forcing frequency. 
The curves present a minimum at the same frequency of 
0.1 Hz, in agreement again with the resonance phenomena. 
However, we observe a difference between the two quanti- 
ties. This underlines that these measured works have not 
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Fig. 3: Distributions of work over one period of Wd (symbols) 
and Wa (solid line) for a forcing at a) 0.05 Hz, b) 0.1 Hz, c) 
0.375 Hz and d) 2 Hz. 



butions of work over a single period. On the figures 3, we 
present the probability density function (PDF) of Wg and 
Wci for different frequencies. One can notice the presence 
of a single bump around the mean value for low frequency 
(Fig. 3a), in agreement with the behavior of the bead. 
The fluctuations reach values, which are larger than six 
times the average injected work. Due to the lack of events 
during an experiment at such low frequency, we are not 
allowed to conclude about the shape of the distribution. 
However, it seems to tend to a gaussian distribution with 
deviations for large fluctuations. Close to the resonance 
frequency (Fig. 3b), the distribution is wide and shifted 
toward a larger value. Then, when the frequency increases 
(Fig. 3c), the distribution shows several peaks. For the 
classical work, the first peak, centered around zero, corre- 
sponds to the work received when the bead does not leave 
its well during the period. The other one corresponds to 
the work done on the system during a single jump of the 
bead. For higher frequencies (Fig. 3d), only the central 
peak remains and increases, the distributions become more 
symmetrical and tends to a gaussian centered on zero since 
the bead hardly leaves its well and thus explores only the 
energy landscape of one well. Although the mean works 
are equal, the distribution of Wg and Wd present signifi- 
cant differences close to the stochastic resonance. In order 
to study in more details these distributions of work and 
heat dissipation, we choose a frequency of external driving 
(0.25 Hz) which ensures a good statistic, by allowing the 
observation of the system over a sufficient number of pe- 
riods, and which, at the same time, produces a non trivial 
distribution (see Fig. 3). We compute the works and the 
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Fig. 4: Evolution of the mean value of works(o, •), dissipation 
(+) and potential variations (x) over an increasing number of 
period n (f = 0.25 Hz). 



dissipation using 1.5 10^ different to on time series which 
spans about 7500 period of the driving. The quantities 
Wsm Wcin and Q„ refer to averages made over n peri- 
ods of the driving. The exchanged heat is defined from 
the first principle of the dynamic: AUq = Wd — Q- The 
figure 4 shows first the linear increase of Wsn, Wdn and 
Q„ with n and second that all these quantities are equals 
because A.Uq vanishes in average (shown by the crosses in 
Fig. 4). The distribution of AUq is plotted in Fig. 5. We 




Fig. 5: Distribution of AUo for n = 1, 2, 4, 
proportional to exp(— Af/o) (/ = 0.25 Hz). 



,12. The lines are 



observe first that these distributions do not depend on the 
number of integration period. This is due to the fact that 
the potential depends only on the position at the end of 
the n-th period, which has the same distribution for all n. 
We notice the classical large exponential wings preceded 
by a sudden decrease in the distribution. This breakdown 
corresponds to the potential energy of the inter-well bar- 
rier and could be understood as follow: The bead can not 
explore so often the high values of the potential as in a 
single well configuration but jumps likely into the other 
well. We can also remark secondary peaks around zero, 
which correspond to a small difference of energy between 
the two wells. Each peak correspond to the probability 



of having a single jump during driving period. For small 
values of AUq, the PFD decreases also exponentially, de- 
scribing the evolution of the bead in a single trap. We 
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Fig. 6: a) Distribution of stochastic work for different number 
of period n = 1, 4, 8 and 12 (/ = 0.25 Hz). Inset: Same data 
in lin-log. b) Normalized symmetry function as function of the 
normalized work for n = 1 (+), 2 (o), 4 (o), 8 (A), 12 (□). 



discuss in the following the PDFs P(X„) of X„ where Xn 
stands for one of the three variables Wgn, Wdn and Qn- 
For these quantities we also study the SSFT which states 
that: 



Sx^ = log( 



PjXn ^ W) , 

P(Xn = -w) ' 



for 



(4) 



where Sx^ is called the symmetry function. It is impor- 
tant to notice that Eq. 4 for Wsn should hold for all n. 
Indeed taking into account that Qn = —At/ -I- Wsn it is 
easy to realize that in this case Wsn/T is just the total en- 
tropy production defined by Seifert [18,23] who has shown 
that that for this quantity the SSFT holds for any inte- 
gration time which is an integer number of periods of the 
forcing. 

We consider first P{Wsn) which is plotted in Fig 6a) 
for various n. The distribution presents a sharp peak for 
n = I, which disappears as soon as n increases. The dis- 
tribution then tends to a gaussian for high values of n. To 
directly test the SSFT, we plot Sws„/ < Wsn > as a fun- 
tion of Wsn/ < Wsn > h^^ Fig. 6b). It is remarkable that 
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straight lines are obtained even for n close to 1, where the 
distribution presents a very complex and unusual shape. 
These curves seem to collapse on the curve y = x. We 
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7: Slope at the origin of the symmetry functions of Wsn 
Wcin (n) and Qn (o) as function of n (/ — 0.25 Hz). 



have looked closer to the slope of Sw^n for small work 
fluctuations. Its evolution is shown on the figure 7. The 
SSFT is not verified at small n but it is rapidly valid for 
increasing number of periods. The convergency turns out 
to be much faster than that observed in a harmonic oscil- 
lator [15] where SSFT is satisfied only for n > 30. This 
is indeed an important point because as we have already 
mentioned Wsn/T corresponds, in this case, to the total 
entropy production, for which SSFT has to be verified for 
all n [23]. Thus the statistical and numerical uncertain- 
ties could explain the fact that in our data and in those 
of ref. [19] the slope of Sw^^ versus Wsn for ti = 1 is not 
exactly one. Indeed in ref. [25] it has been shown that 
even in the gaussian case it is required a very large sta- 
tistical accuracy in the calculation of the total entropy to 
satisfy SSFT for small n. In our case this accuracy has 
to be larger than in the gaussian case because of the ex- 
tremely complex shape of the P{Wsn) for n=l. For the 
classical definition of the work Wdn, the fluctuations are 
larger than those of Wsn (Fig. 8a). Again, the distribu- 
tions tend to a gaussian for large n (inset of Fig. 8a). On 
Fig. 8b), we have plotted the normalized symmetry func- 
tion of Wcin- We can see that the curves are close to the 
line of slope one. For high values of work, the dispersion 
of the data increases due to the lack of events. The slopes 
at the origin are shown in Fig. 7. We can see an increase 
of the slope near zero and then an oscillation. Although 
we can not average over more periods, the slope seems to 
tends toward 1 as expected by the SSFT. The strong anal- 
ogy existing between the convergency of Wsn and Wd is 
probably due to the fact that, as already mentioned, they 
differ only for a boundary term that rapidly goes to zero 
when n is increased. Now, we consider the heat fluctu- 
ations. The PDFs show even larger fluctuations (Fig 9). 
We can notice flrst that the PDFs decrease exponentially 
in the asymptotic limit. This is due to the exponential 
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Fig. 8: a) Distribution of classical work Wd for different num- 
bers of period n = 1, 2, 4, 8 and 12 (/ = 0.25 Hz). Inset: 
Same data in lin-log. b) Symmetry function as function of the 
normalized work (same symbols as in Fig. 6.) 



tails of the distribution of AUq that become predominant. 
This behavior is less pronounced for higher values of n. 
Near the maximum, these distributions tend to be gaus- 
sian as shown by the exponential fit (dashed line) . On the 
figure 9b), the symmetry function of the dissipated heat 
is shown as function of the dissipated heat. We can define 
two regions, in the first one, the curve increases linearly 
and then bents. On the inset of Fig. 9b), which shows 
the normalized symmetry function of Qn, the two differ- 
ent regimes appear clearly. First, the slope of the linear 
part increases with n and goes asymptotically toward 1, 
as shown on Fig. 7. Only the linear part for Qn < 3 has 
been fitted. Starting near 0.6 for one period, it grows to- 
ward 1 more slowly than for the classical work. This shows 
that the SSFT holds asymptotically when n increases for 
Qn ^ {Qn)- On the other regime, one could expect the 
curve to reach the value 2 from the shape of the expo- 
nential tails. However, their size are too small to create 
a horizontal asymptote. Still, when normalizing the data 
by the average of Qn (inset of Fig. 9b), the data begin to 
explore the second regime around a value compatible with 
the value 2, that would need longer measurements (over 
more period) to be confirmed. 
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Fig. 9: a) Distribution of the dissipation for different numbers 
of period n = 1, 2, 4, 8 and 12. The straight line is proportional 
to exp{—Q,T). The dashed line is a gaussian fit of the PDF for 
n = 12. Inset: Same data in lin-Iog . b) Symmetry function 
of the dissipation for n = 1, 2, 4, 8 and 12. Inset: normalized 
symmetry function (/ = 0.25 Hz), (same symbols as in Fig. 
6.) 



In conclusion, we have experimentally investigated the 
power injected in a bistable colloidal system by an ex- 
ternal oscillating force. We find that the injected power 
in the stochastic resonance regime presents a maximum 
when the frequency of the driving force corresponds to 
half of the Kramers' rate. We have compared the stochas- 
tic work and the classical work for different frequencies and 
shown that although the average values of each work are 
equal, the distributions of work reveal some differences. 
They present large tail toward negative values and their 
shape are non-gaussian close to the resonance when aver- 
aging over a single period. Analyzing in more details the 
distribution of work and dissipated heat, we have tested 
the validity of SSFT for a non-linear potential. We have 
shown that FT rapidly converge to the asymptotic value 
for rather small n. The fact that for the total entropy 
WsnlT^ SSFT is not satisfied exactly for small n is cer- 
tainly due to statistical and numerical inaccuracy. We 



have also shown that the SSFT is only valid for small val- 
ues of dissipated heat compared to the mean value at long 
time. It would be interesting to change the symmetry of 
the driving force cycle [26] to explore the limit of the FT. 
We acknowledge useful discussion with R. Benzi, K. 
Gawcdzki and S. Joubaud. This work has been partially 
supported by ANR-05-BLAN-0 105-01. 
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